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1. Introduction 

There has been considerable interest in understanding the electronic properties of a 
single layer of graphene ever since its experimental realisation. Experimental and 
theoretical studies have shown that the nature of the quasi-particles in these two- 
dimensional systems is very different from that of the standard two-dimensional electron 
gas (2DEG) which has been extensively studied. Graphene has a honeycomb lattice 
of carbon atoms. The quasi-particles in graphene have a band structure in which 
electron and hole bands touch at two points in the Brillouin zone. At these Dirac 
points the quasi-particles obey the massless Dirac equation: they behave as massless 
particles with a linear dispersion relation = vk (with the characteristic velocity 
v ~ 10 6 ms _1 ). This behaviour gives rise to a host of new and unusual phenomena 
such as anomalous quantum Hall effects and a n Berry phase [HE]- This 2D Dirac-like 
spectrum has been confirmed by measurements of de Haas-van Alphen and Shubnikov- 
de Haas (SdH) oscillations [3], where magnetic oscillations appear due to the interplay 
of the Landau levels with the Fermi energy, and are important tools in the investigation 
of the Fermi surface and electron transport. In a standard 2DEG an artificially 
created periodic potential in the sub-micron range leads to the appearance of Weiss 
oscillations in the magnetoresistance. Such electrical modulation of the 2D system can 
be achieved by depositing an array of parallel metallic strips on the surface or through 
two interfering laser beams [H EJ [6]. Weiss oscillations can be explained in terms of 
the commensurability of the electron cyclotron diameter at the Fermi energy and the 
period of the electric modulation. These oscillations were found to be periodic in the 
inverse magnetic field |5j El [7] . It is therefore interesting to study the affect on the Dirac 
electrons of electrical modulation of a graphene layer. In this work we study the effects 
of modulation on the bandwidth (A) and the density of states (DOS) of the Dirac 
electrons in graphene. These quantities are essential prerequisites for understanding 
properties such as electron transport, thermodynamic behaviour etc. 

In section [2], we present the formulation of the problem. Section [3] contains the 
calculation of the density of states whereas in section [4] we discuss the bandwidth for 
electrically modulated graphene including an asymptotic and classical description. The 
conclusions are in section [5j 

2. Formulation 

We consider two-dimensional Dirac electrons in graphene moving in the x-y-plane. 
The magnetic field, B = (0,0, B), is applied along the z-direction perpendicular to 
the graphene plane. This system is subjected to a weak electric modulation along the 
x-direction. Using the Landau gauge we write the vector potential as A = (0,Bx,0). 
The two-dimensional Dirac like Hamiltonian for a single electron in the Landau gauge 
is (using /i = c=l)[D[2j[8] 




(1) 
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where er = {o~ x ,o~ y } are the Pauli matrices and v is the magnitude of the electron 
velocity. The complete Hamiltonian of our system may be written as 

H = H + U{x) , (2) 

where Hq is the unmodulated Hamiltonian and U(x) represents the periodic modulation 
along the x-direction modelled as 

U(x) = V Q cos(Kx) , (3) 

where K = 2ir/a, and a and Vq are the period and amplitude of the modulation 
respectively. Without modulation the Landau level energies are given by 

e(n) = uj g y/n, (4) 



where u g = vyleB is the cyclotron frequency of the graphene electrons and n is an 
integer. Note that the Landau level spectrum for Dirac electrons is significantly different 
from that in a conventional 2DEG where e{n) = oj c {n + |) and uj c = eB/m is the 
cyclotron frequency. 

The eigenfunctions without modulation are given by [8] 

* k(r ,J^(-^[(x + x )/i]\ 



$n0) = / ——— H n (x) , (6) 



where 



I = (eB)- 1 / 2 is the magnetic length, the ^/-dimension of the graphene 

layer and H n (x) are the Hermite polynomials. 

As we are considering weak modulation such that Vo is smaller than the Landau 
level separation we can apply standard perturbation theory to determine the first order 
correction to the unmodulated energy eigenvalues 

/oo r-Ly 
dx / dy^ ky (r)U(x)^ n , ky (r) (7) 
-oo JO 

with the result [8] 

AE n , ky = \V Q cos(ivx )e-"/ 2 [L n ( M ) + L n ^{u)} (8) 

where u = K 2 £ 2 /2 and L n (u) are Laguerre polynomials. Hence the energy eigenvalues 
in the presence of modulation are 

e(n, x Q ) = e(n) + AE n<ky = u g y/n + \F n \ cos(Kx Q ) (9) 

with \F n \ = ^V e~ u ^ 2 [L n (u) + L n -i(u)]. We observe that the degeneracy of the Landau 
level spectrum of the unmodulated system with respect to k y (and x ) is lifted in the 
presence of modulation. The formerly sharp Landau levels broaden into bands whose 
widths, ~ |F n |, oscillate as a function of n since L n (u) are oscillatory functions of the 
index n. At this stage we can compare the energy spectrum of Dirac electrons with that 
of standard electrons in the same system. The differences are: 
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• the standard electron unperturbed energy eigenvalues depend linearly on both the 
magnetic field and the quantum number n whereas for Dirac electrons they depends 
on the square root of both. 

• in graphene we have the average of two successive Laguerre polynomials [L n (u) + 
L n -i(u)]/2 while for standard electrons we have a single term, L n (u). 

These differences will give different results for the density of states and the band 
width, as we show in the next section. Note that for the weak electric modulation 
case under consideration the quantum numbers n can be referred to as the magnetic 
Landau band indices and are equivalent to the Landau level quantum number n for 
the unmodulated system. In the presence of weak electric modulation, the band width 
of the magnetic Landau bands depends on the index n. Thus the electric modulation 
induced broadening of the energy spectrum is non-uniform, a feature which will be of 
significance in understanding the behaviour of Dirac electrons in modulated graphene. 



3. The Density of States (DOS) 

It is well known that in the absence of modulation the DOS consists of a series of delta 
functions at energies equal to e(n). The addition of a weak spatially periodic electric 
modulation, however, modifies the formerly delta function like DOS by broadening the 
singularities at the energies (e(n)) into bands. The density of states is given by 

D(e) = jJ2 6 ( £ - £ ^ y ) , (10) 

nky 

where the sum on n extends over all occupied Landau levels and A is the area of then 
sample. By using the energy eigenvalues given in we can express D(e) as: 

D ( £ ) = 2 ^pJ2f dx 5(e-e n -\F n \ C os(Kx )), (11) 

where e n = uj^^/n, and a factor 2 is due to spin degeneracy. Evaluation of the x — integral 
in the above equation yields the zero temperature density of states of the modulated 
two-dimensional Dirac electrons: 

D ( e ) = -m E , 1 = e d F "l ~\ £ - £ n\), (12) 



where Q(x) is the Heaviside unit step function. Here we can see that the one-dimensional 
van Hove singularities of the inverse square-root type appear at the low and high energy 
edges of the broadened Landau bands, forming a double peak like structure. 

The zero temperature DOS given by (I12p is shown graphically in figure [J as a 
function of energy, using the following parameters [8]: v ~ 10 6 ms" 1 , hd = 3x 10 15 m~ 2 , 
a = 350 nm, Vq = 0.35 meV, and kp = (27mD) 1//2 being the Fermi wave number 
of the unmodulated system in the absence of a magnetic field. The origin of both 
Weiss and of Shubnikov-de Haas (SdH) oscillations is immediately apparent. The short 
period, high amplitude oscillation is the Landau level structure which gives rise to SdH 
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Figure 1. The dimcnsionless density of states, D/Dq, in a periodically modulated 
graphene as a function of energy (a & b have different energy ranges only, all other 
parameters are the same) for fixed value of magnetic field B = 0.35 T. 



oscillations whereas the apparent darker longer period oscillation (actually the minima 
in the DoS for each Landau level) is associated with the Weiss oscillations. This is 
largely a consequence of the oscillatory factor |-Fjv|, which has been shown to exhibit 
commensurability oscillations. Our basic density of states spectrum is exactly the same 
as shown by others [9] without modulation, as is the Weiss period, albeit with an extra 
modulation. 
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4. The Bandwidth (A) 

To better appreciate the modulation of the amplitude of Weiss oscillations we plot the 
bandwidth as a function of the magnetic field in figure [21 The width of the n th Landau 
level is given as 

A = 2|F Jv i = V exp-3|L n («) + L n _ 1 (w)| . (13) 




B(T) 

Figure 2. Bandwidth due to periodic electric modulation 2D graphene as a function 
of magnetic field. The dotted line is the asymptotic behaviour and solid line represents 
the exact behaviour of the width. 

This is clearly different from the standard electron result jH [5j [6]. The bandwidth is 
plotted for n = n-p where n-p = Ep/u 2 is the Landau level index at the Fermi energy. 
Contrast this with rip = E-p/oj c — | and to c = eB/m for a standard 2DEG. For the 
low magnetic fields under consideration, the graphene results for Dirac electrons are the 
same in phase and amplitude as those for standard electrons. Moreover, we have found 
the maxima and minima of the bandwidth at the same points as for the case of standard 
electrons jU [7] . 

Asymptotic Expression 

An asymptotic expression for the bandwidth can be obtained by using the following 
expression for the Laguerre polynomials in the limit of large n as 

e~ u/2 L n (u) -> (tt W) ~ 1/4 cos (2(nu) 1/2 - f ) . (14) 

Substituting the asymptotic expression given by (fl~4"|) into (1T31 yields the asymptotic 
expression for bandwidth 

A = V (7rW)~ 1/4 cos {\{u/n) l l 2 ) cos (2(nu) 1/2 - f) (15) 
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where we have rewritten (fT5|) containing w = K 2 £ 2 /2 in terms of the ratio of the semi- 
classical orbital radius R g and the modulation period a. This expression can be easily 
understood by analogy with the beating of 2 oscillators of similar frequencies: the first 
cosine term is the amplitude of the beat. This extra modulation of the bandwidth is 
the most significant difference between the DOS of a normal 2DEG and graphene. Note 
that for large n it approaches unity as {u/n)~ l . 

For large n the level spacing goes as uj g (n 1 / 2 — [n — l) 1 ^ 2 ) — > uj g \n~ 1 / 2 and the 
width goes as (ir 2 nu) 1//4 , apart from the modulation. There is therefore a value of n at 
which the width becomes equal to the spacing and the perturbation theory is no longer 
valid. This occurs when 

vr 4 

1 -2„,, ,4 _ " „,4 



"max = Ye * VW g = T^yl^ eB . (17) 

For a fixed electron density this suggests a minimum value for the magnetic field B 
below which it is necessary to carry out a more sophisticated analysis. Note that this 
argument also applies to any other calculation [8J which treats the modulation as a 
perturbation. 

4-2. Classical description 

We now give a classical explanation of the asymptotic expression of bandwidth obtained 
in ffl6|) which is essentially a large n expression. The classical equations of motion along 
the x and y directions are 

x(t) = x Q + R g sin(u; g t + ip) (18a) 
y{t) = Vo + Rg cos(u g t + cp) (18°) 
respectively, where R g is the radius of the orbit, xo and yo are the centre coordinates 
and (p is phase factor. Note that this approach is valid for both graphene and parabolic 
systems apart from the fact that the orbital radius R g scales as Ep in the graphene 

1/2 

case but as E F otherwise. Without loss of generality we may take ip = 0. Thus the 
increase in the average energy of the cyclotron motion due to the electric modulation is 
evaluated as 

i r+to/2 

AE(x ) = - V cos(Kx(t))dt (19) 

*0 J -to/2 

where to is the period of the orbit. This result is valid to the same order as (J7|). 
Substituting x(t) yields 

AE(x ) = V J (KR g ) cos(Kx ) (20) 

with Jq(z) the Bessel function of zero order. For 2irR g > a, one can replace the Bessel 
function Jq by a cosine function as 
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with the result 

which is almost the same as obtained in ffl"6]) in the limit of large n. The significant 
difference between ( Ti~6l) and ( 1221) is the extra cosine factor which has its origin in the 
mixing between consecutive Landau levels, or Laguerre functions, which is something 
with no classical equivalent. This reveals itself in ( fl6l) in that it has not been possible 
to rewrite (TT5I) in terms of purely classical quantities: there is a residual n in the extra 
factor. 

5. Conclusions 

In this work we have analysed the band spectrum of graphene with a magnetic field 
perpendicular to the graphene layer and a unidirectional electric modulation. We have 
determined the density of electronic states and the bandwidth of each level. We have also 
considered the asymptotic expression for the bandwidth and its relation to a classical 
description and have noted a quantum correction to the classical behaviour. To highlight 
the effects of modulation on the density of states and bandwidth, we have plotted these 
quantities for experimentally relevant parameters. 
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